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ABSTRACT 

An approximate mathematical analysis of supersonic inlet dynamics is presented. 
The subsonic duct is subdivided into constant -area sections which are represented by 
one -dimensional wave equations. The movable normal shock is used as the upstream 
boundary condition; and a choked station is assumed lor the downstream boundary con- 
dition. The analysis leads to a closed “form matrix solution for the shock -posit ion 
transfer function* Analytical results are presented and compared with experimental 
data. 
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A MATHEMATICAL ANALYSIS OF SUPERSONIC INLET DYNAMICS 

by Ross G. Wiiloh 
Lewis Research Center 


SUMMARY 

An approximate mathematical analysis of supersonic inlet dynamics is presented. 
The subsonic portion of the inlet is represented by one-dimensional wave equations which 
are integrable along constant -area portions of the subsonic duct. The complete subsonic 
duct is thus represented by subdividing it into constant -area sections. The movable nor- 
mal shock is used as the upstream boundary condition, and a choked station is assumed 
for the downstream boundary condition. The analysis leads to closed-form matrix solu- 
tions for shock position and pressure frequency response. Results from the analysis ax'e 
presented and compared with test data from an experimental inlet. 


INTRODUCTION 

Conventional propulsion systems for supersonic aircraft consist of either a turbojet 
or turbofan engine and a supersonic inlet. The inlet is intended to decelerate the super- 
sonic air stream to velocities suitable for the engine and, in so doing, convert a portion 
of the free-stream kinetic energy to static pressure. In a mixed-compression inlet, 
supersonic deceleration is accomplished in a converging section of the inlet duct and sub- 
sonic deceleration in a diverging section. The two regimes are separated by a normal 
shock just downstream of the minimum -area portion of the duct. 

Inlet pressure recovery is highest when the normal shock is positioned at, or just 
downstream of, the minimum -area point, if, however, the shockwave moves upstream 
of the minimum -area position, it generally becomes unstable and continues to move until 
it stands in front of the inlet. When the shock is in front of the inlet, engine weight flow 
and thrust are sharply reduced, and vehicle drag is increased. 

Inlets with variable geometry are used to achieve high performance for aircraft 
operating over a range of supersonic speeds. The use of variable inlet geometry does 
not in itself guarantee the desired performance, since variable geometry is an asset only 
when it is properly positioned. Control systems are needed to vary the inlet geometry 



not only for steady-state flight requirements, but also to minimize the effects of system 
disturbances arising from gusts, aircraft maneuvers, and changes in engine operation. 
The design of a high-performance, shock-position-control system requires frequency 
response data or a mathematical description of shock-position dynamics. The purpose 
of this report is to develop an approximate mathematical analysis for the dynamics of 
supersonic inlets. A means for obtaining shock-position frequency response and a 
mathematical representation suitable for analog simulation are presented. 

The theory of shock motion in channel flow is examined in a number of ways in the 
available literature. One of the earlier reports in the area (ref. 1) is concerned with 
the stability of shock waves in channel flow. In reference 2, the concept of a Helmholtz 
resonator is applied to an inlet. In reference 3, the method of coordinate perturbation is 
applied to the interaction of a disturbance and a shock wave in diffuser flow. In refer- 
ence 4, a Helmholtz resonator is combined with a single capacitive lump. The analysis 
of this reference is carried out in detail, and the results are compared with experimen- 
tal data. In addition, a discussion of the application of the method of characteristics is 
presented. In reference 5, the effects of pressure disturbances just downstream of the 
normal shock on shock position are investigated in a linearized analysis. 

The experimental investigation of inlet dynamics presented in reference 6 indicates 
the presence of multiple resonances typical of distributed -parameter systems. Any 
analysis which predicts such multiple resonances must contain a duct representation 
which includes distributed effects. The number of lumps needed with conventional lump- 
ing techniques rises rapidly with any requirement for accuracy beyond the first reso- 
nance. The method presented in this report combines a set of linearized equations 
across the normal shock with an exact solution of the linearized wave equation. This 
treatment of the subsonic duct avoids the complexity of the method of characteristics 
while still predicting the multiple resonances obtained in reference 6. The analysis pre- 
sented is more exact than conventional lumped -parameter techniques, and for many 
problems, is no more complicated in application. 


MATHEMATICAL ANALYSIS 

An idealized one -dimensional inlet is shown schematically in figure 1. The sym- 
bols used in this figure and throughout the text are defined in appendix A. A linearized 
dynamic analysis relating shock position to adjacent parameters is presented in appen- 
dix B. Linearized wave equations for the subsonic portion of the duct are developed in 
appendix C. The shock-position equations become the upstream boundary conditions for 
the wave equations. A choked exit station was assumed for the downstream boundary 
condition. The mathematical formulations are based on the selection of flow, total 
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Figure L - Idealized mixed-compression inlet. 


pressure, and entropy as the independent variables to describe system behavior. Al- 
though other parameters may be used, those selected are constant in the steady state, 
since friction losses are neglected, and thus considerably simplify the linearized analy- 
sis. 


Shock- Position Dynamics 

The equations for shock position are developed in appendix B. Continuity, momen- 
tum, and energy balances are taken across the normal shock and modified to account for 
a moving shock. With the use of the assumptions of perfect -gas relations, a constant 
specific -heat ratio of 1. 4, and negligible change in duct area across the normal shock, 
the equations are Linearized in terms of the selected variables. It is shown in appen- 
dix B that 
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where the subscripts 1 and 2 designate fixed stations just upstream and just down- 
stream of the shock. 

The dynamic relations for shock position given in equations (1) to (3) are valid for 
upstream and downstream perturbations. If only downstream disturbances are consid- 
ered, AP t j, AW j, and ASj are zero; and equations (1) to (3) can be written in terms 
of the Laplace transform as 


where 
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Subsonic-Duct Equations 

The analysis of the subsonic duct is based on a linearization of the compressible - 
flow equations written in terms of weight flow, total pressure, and entropy. Using the 
assumptions of perfect-gas relations and an inviscid non -heat -conducting fluid, the 
following set of equations are developed in appendix C as equations (95): 
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After linearization, if area changes are assumed to be slow, these expressions reduce 
to equations (99) of appendix C. 
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The equations of the above set are in the form of wave equations, and their solutions 
may be written in terms of a time shift, a is thus shown in appendix C that the fluid 
properties between two stations, 1 and 2, of a subsonic duct can be related in Laplace 
transform notation by 


P 2 + otir 2 + <3S 2 = e‘ as (p 1 + a&j + CO 

P 2 - aW 2 - ySg = e TS (P 1 - aW x - ySj) (8) 

S 2 = e" 08 ^) (9) 
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where cr = L/{v + a), t * L/(a - v), and 6 = L/v. The integration over the coordinate 
x, used in arriving at equations (7) to (9) is based on a constant -area duct. These equa- 
tions must thus be regarded as a solution for the case of a constant -area duct. 

Equation (7) represents the propagation of a wave composed of a linear combination 
of total pressure, weight flow, and entropy. The wave travels downstream from sta- 
tion 1 to station 2 at a velocity equal to the sum of the fluid velocity and the speed of 
sound. Equation (8) represents the propagation of a wave composed of another linear 
combination of the same variables. This wave travels upstream from the second station 
to the first at a velocity equal to the speed of sound less the flow velocity. Equation (9) 
represents an entropy wave traveling with the fluid at the fluid velocity. 

Since the integration necessary to obtain equations (7) to (9) is valid only for a 
constant -area duct, variations in duct area are included in the analysis by dividing the 
duct into sections. The Mach number in each section is averaged by summing one-half 
of the initial and final values along the section. The parameters a, j3, and y are then 
computed from the averaged Mach number. The computation of the corresponding delay 
times a, t, and 6 is also based on average properties in the section. With the evalu- 
ation of these constants, equations (7) to (9) form a set which describes each subsonic 
section. 

Equations (4) to (9) account for shock and sub sonic -duct dynamics. The set thus 
provides a basic mathematical model which can be applied to a supersonic inlet. For the 
solution of a specific problem, however, the disturbance must be described in terms of 
the dependent variables, and an equation for the termination of the final subsonic section 
must be provided. 


Disturbance Equations 

Although the basic mathematical formulation presented will accommodate various 
disturbances, consideration will be restricted to flow perturbations introduced in the 
subsonic duct. Figure 2 illustrates such a disturbance and the notation associated with 
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Figure 2. - Flow disturbance In subsonic duct. 
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it. If the section preceding the disturbance is terminated at station n just before the 
disturbance, the following equations relating fluid conditions before the disturbance to 
those just downstream can be written in Laplace transform notation as 


rV k'nJ 

P n = P n+1 

(10) 

■ *n + l + 

(ID 

rsv 

S = S t 
n n+ 1 

(12) 


Exit Boundary Condition 

The equations applied at the end of the last section depend on the problem boundary 
conditions. E it is assumed that the inlet is terminated with a choked orifice, the flow 
can be described by 
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(13) 


where M = 1. 0 for a choked orifice, and it is assumed that k = 1.4. Linearizing equa- 
tion (13) and simplifying with 
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Equation Summary 


The shock-position or upstream boundary equations (4) to (6), the subsonic duct 
equations (7) to (9), the disturbance equations (10) to (12), and the exit boundary condi- 
tion, equation (14), form a complete set describing shock position to perturbations in 
duct flow. Rewriting equation (8) as 

?1 - aftj - ySj = e" TS (p 2 - aW 2 - yS 2 ) (15) 

allows the direct simulation of the set of equations. In this form, the equations may be 
represented on a digital computer with one of the simulation languages such as MIMIC 
(ref. 7) 


MATRIX SOLUTION OF EQUATIONS 
Mathematical Formulation 


Even though direct simulation of the complete set of equations is difficult because of 
the algebraic complexity, a closed-form matrix solution for the shock-position transfer 
function can be obtained. This solution can be programmed for the numerical evaluation 
of the desired transfer function. 

In matrix form, the shock-position expressions, equations (4) to (6), become 
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t 2 = gx 


(17) 


If station 2, just downstream of the shock, is separated from station 3 by a section of 
subsonic duct, equations (7) to (9) can be used to relate conditions at the two stations, 
hi matrix notation these equations become 
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or 


B 2 T 3 " D 2 B 2 T 2 


(19) 


]ii a similar manner, stations k and k + 1 are related by 


B k T k + l =D k B k T k 


( 20 ) 


In matrix form the expressions describing the disturbance flow, equations (10) to (12), 
become 
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( 21 ) 


or 
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( 22 ) 


while equation (14) for the exit boundary condition becomes 


0 

0 


0 

0 



!P J 


fV 

", 


M 


= 0 


(23) 


or 


FT- = 0 


(24) 
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Equations (17), (19), (20), (22), and (24) can be combined to represent the inlet dynamics. 


t 2 = gx 
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With these definitions, equations (25) can be written as 


T n ' E n-I E n-2 ' ' • E 3 E 2 GX 
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If the products of E matrices are replaced by 


K = E n-l E n-2 ‘ ‘ ‘ E 3 E 2 
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and 


N = E i-l E i-2- 


E n+2 E n+l 


equations (26) become 
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Equations (27) can be combined into 


FNKGX = FNT t 


(27) 


(28) 


Since FNKG has only a single nonzero element because F is a single-row matrix and 
G is a single -column matrix, equation (28) can be written as 


F l,a N a,/3 K j3,y G y, 1 X 1 ~ F l, a N a, 2 T D, 2 


(29) 


where the summation convention (summation across repeated subscripts) is assumed. 
As =X g and 2 = W^, then 


F i,Ai 


„N, 


3 K 0 G 


W D 1, o' a,Fj3, y Y, 1 


(30) 


The elements of the N, K, and G matrices are functions of s. Equation (30) thus rep- 
resents the transfer function for the response of shock position to the flow disturbance 
assumed. For s - ioj, the frequency response, in amplitude and phase, is easily and 
quickly calculated numerically by making the operations indicated by equation (30) in 
complex algebra. Substituting these results for shock-position response into the first of 
equations (25) gives the response for the elements of T 2 ; substituting these results into 
the second of equations (25) gives the response for the elements of T^, etc. Thus, the 
frequency response for perturbations in total pressure, flow rate, and entropy are ob- 
tained at all stations. 
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The perturbations in other inlet variables, such as static pressure, velocity, etc. , 
can be obtained at each station from those in total pressure, flow rate, and entropy. An 
example of this is given later (eq. (34)). 


Comparison of Analysis with Test Data 

Frequency response tests were conducted on a 48 -centimeter mixed- compression 
inlet in the Lewis Research Center 10- by 10 -foot supersonic wind tunnel. A detailed 
description of the inlet is reported in reference 8, and the data resulting from the fre- 
quency response test program and a description of the disturbance device are reported 
in reference 6. To verify the analytical approach presented in this report, the techniques 
have been applied to this 48 -centimeter inlet. 



Figure 3. - Mixed-compression, 48 -centimeter iniet, 


A schematic diagram of the 48 -centimeter mixed-compression inlet is presented in 
figure 3. The variable-area portion of the inlet, beginning at the shock and extending to 
the disturbance generator, was divided into three sections of equal length. The re- 
mainder of the inlet, extending from the disturbance generator to the choked exit, was 
represented by a single section. With this sectioning, equations (26) become 


T 6 . E 4 E 3 E 2 GX 


T e - T s - t d 


t 7 = e 6 t 6 


FT, - 0 


( 31 ) 
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or, after appropriate manipulation, 


FEgE 4 E 3 E 2 GX = FEgT D (32) 

A schematic diagram of the inlet is presented in figure 4 to indicate the sectioning tech- 
niques and the location of significant problem variables. Table I contains numerical 
values associated with the tests of the 48 -centimeter inlet. In addition, coefficients, 
based on steady- state, one -dimensional flow calculations averaged across each section, 
are presented for use in the evaluation of equations {31) and (32). 


;i Choked exit 



Figure 4. - Schematic indicating station numbers and associated parameter subscripts. 


TABLE L - NUMERICAL VALUES FOR 48 -CENTIME TER -INLET 

COEFFICIENTS BETWEEN STATIONS 

jcoefficients for shock conditions: = 0. 3133; C 2 = 5. 483x10 sec; 

C 3 = 2. 013xl0' 3 sec. Test conditions: total pressure, 2534.4 1b/ft 2 
(121 130 N/m 2 ); total temperature, 689.3° R (382.9 K); area ratio 
A’/A, 0.4165 ft" 1 (1.366 m 1 ).] 


Station 

Coefficient 

Transportation time, 

sec 


a 

JS 

y 

0 

6 

T 

2 to 3 

0.9225 

0. 3597 

|0. 1675 

0. 2961X10' 3 

0. 7034xl0~ 3 

l,872xMT 3 

3 to 4 

* 6752 

,2418 

* 1440 

, 3309 

* 9845 

1,009 

4 to 5 

. 4956 

. 1673 

* 1159 

,4573 

1* 716 

.9789 

6 to 7 

.4222 

,1392 

- 1021 

4,075 

17* 44 

7.690 
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A digital program to perform the necessary operations was written in FORTRAN IV. 
The coefficients required for the various matrices are placed in storage and subsequently 
used to calculate the matrix elements. The matrices are formed as arrays from these 
elements, and the matrix multiplications of equation (32) are then performed. Both the 
G matrix and the E matrices contain the Laplace transform operator s which is re- 
placed by io> for numerical evaluation of the frequency response. These matrices and 
their products contain complex numbers, and the program is thus written in complex 
algebra. 

The general scheme is to select a frequency value and compute the amplitude 
and phase of the complex number 


X s< iw) , , 

= H < iw > = *>3 + id 3 

W D (io>) 


(33) 


which results from the evaluation of equation (32). The frequency value is then incre- 
mented, and the procedure is repeated. The resulting tabulation yields data for a plot 
of the shock-position -to -disturbance -flow frequency response. 

After evaluation, the shock-position-to-disturbance-flow transfer function was used 
with equations (31) to evaluate the response of the pressure immediately downstream of 
the shock. Equations (31), however, are written in terms of total pressure, while avail- 
able test data is in static pressure. The equation relating static pressure to the system 
dependent variables at some station j is 

M 2 + 5 AP T 7M 2 AW . M 2 AS ^ 34 j 

s(l -M 2 ) P T 5(l-M 2 )w 5(l -M 2 ) R 



Equation (34) was used with equations (31) in the actual evaluation of the static -pressure 
responses presented. 

Figures 5 and 6 contain the results of the computations plotted with the correspond- 
ing test data from reference 6. The amplitude ratios are plotted in normalized form; 
that is, the amplitude ratio measured at the test frequency is divided by the amplitude 
ratio determined at a low reference frequency. 

Figure 5 contains the shoek-position-to-disturbanee-flow data. The figure indicates 
the good agreement that was achieved between the experimental and analytical shock- 
position results. 

Figure 6 contains response plots for static pressure just past the normal shock to 
disturbance flow. Again, both test data and analytical results are presented. Although 
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the test results have slightly less peaking than predicted by the analysis, agreement is 
reasonably good. 


Model Simplification 


As previously indicated, the complexity of the waves arriving at the nodes between 
the sections of the subsonic duct complicates the simulation of the system equations. An 
examination of the subsonic -duct equations (7) to (9) reveals that the model complexity 
would be greatly reduced if constant entropy were assumed. Although such an assump- 
tion is in no way required for the evaluation of inlet transfer functions through numerical 
calculations, it is required for any reasonably simple simulation of the inlet. 

E constant entropy is assumed, a change in the zero -frequency shock-position gain 
results. This gain change arises from the system boundary -condition equations. Apply- 
ing the final-value theorem to equations (4) and (6) allows evaluation of the steady -state- 
pressure-to-shoek-position and entropy -to -shock -position zero-frequency gains. 


X, 
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The exit boundary condition is given in equation (14) as 



This equation, together with the zero -frequency gains at the inlet boundary, implies that 


W 


E 


= 1 


w— 0 


E, however, S is set equal to zero throughout the problem, equation (14) becomes 
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which implies that 



W 


E 


w-»0 


1 

6 


Thus, the assumption of AS = 0 results in a system gain change. To avoid this change, 
the simplified model assumed that equation (14) could be modified to 


which is equivalent to assuming constant temperature at the exit boundary and constant 
entropy elsewhere in the system. 

To evaluate the effects of the model simplification, the frequency response of the 
four -section model was numerically evaluated with and without the simplifications. Fig- 
ure 7 contains plots of the numerical data. Examination of the curves plotted in figure 7 


S’ 


i Two -sect ion model with 

constant entropy 

— Four-section model with 
constant entropy 




Figure 7, - Shock -posit ion response to flow perturbations. Com- 
puted from mathematical model with and without simplications. 
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indicates that for the inlet considered herein, the simplified model yields good agreement 
in the frequency range of interest. 

To investigate the possibility of further simplification, the frequency response of the 
simplified system was evaluated with a two-section representation where the variable- 
area portion of the duct was represented by a single section. The results of these calcu- 
lations are also plotted in figure 7. The additional simplification did not significantly 
alter the results. 


Analog Simulation 

The matrix methods presented yield frequency response data for the inlet dynamics. 
To secure results in the time domain and to provide a working model for control-system 
synthesis, an analog simulation was investigated. Since the numerical evaluation of the 
system equations indicated that a two -section, constant -entropy representation gives 
reasonable agreement with the test data, this representation was made the basis for an 
analog simulation. The system transportation times (<j and r) were represented with 
the fourth-order Pade approximations of reference 9. Frequency response data obtained 


Analog simulation with 
constant entropy 



Figure & - Shock -position response to flow pertu rations. Com- 
puted from analog simulation and from complete equation set 
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from the analog simulation are plotted in figure 8 along with the digital solution for the 
complete set of equations. The agreement shown in this figure is excellent to frequen- 
cies as high as 200 hertz. The growing discrepancies beyond 200 hertz can be attributed 
to the Fade approximations whose accuracy deteriorates with rising frequency. 

Figure 9 contains a shock -position, time -history response to a step change in bleed 
flow, obtained from an analog simulation with Pade approximations used for the delay- 
time representation. Once established, the analog simulation can be used as a design 
tool for inlet shock-position-control design and experimentation. 


20xl0 3 



Time, t, sec 

Figure 9. - Shock-position response to step change in bleed flow, 
analog computer simulation* 

CONCLUDING REMARKS 

A theory for supersonic inlet dynamics has been presented based primarily on the 
closed-form solution of the linearized one -dimensional distributed -parameter wave 
equations. For transfer functions, closed -form expressions were obtained in terms of 
matrix operations (no matrix inversion required). For numerical evaluation of frequency 
responses, these matrix operations can be simply programmed on a digital computer in 
complex algebra. The results of the theory were found to be in excellent agreement with 
data from a 48 -centimeter inlet, Various further simplifications to the theory were 
evaluated. The theory developed herein can be applied also to analog computer simula- 
tion. 
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In the application of the methods presented to the prediction of shock-position dy- 
namics, difficulty arises in the selection of an appropriate value of A'/A (i. e. , 

(1/A) (dA/dx) in the vicinity of the normal shock) to be used. Because of boundary layer 
effects and flow distortion in the throat region, shock dynamics are not uniquely depen- 
dent on the geometric A f /A. The value used in the analysis should thus include the in- 
fluence of the boundary layer as well as the inlet geometry. Further work is required to 
allow the analytical determination of an effective A'/A parameter. 

In general, results obtained from both the full and simplified models compare well 
with test data obtained from experimental inlets. Further work must be done to deter- 
mine the regions to which the model can be successfully applied, and to include the ef- 
fects of performance bleed and other variations on the basic inlet. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 28, 1968, 

720-03-01-04-22. 
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APPENDIX A 


SYMBOLS 


A 

2 2 

cross-sectional area, ft ; m 

K 

matrix 

A' 

rate of change of cross-sectional 

k 

ratio of specific heats 


p 0 

area along duct, ft /ft; m /m 

L 

length, ft; m 

a 

speed of sound, ft/sec; m/sec 

L>T 

duct length from normal shock to 

B 

coefficient matrix 


exit station, ft; m 

b n 

constant coefficients {n = integer). 

M 

Mach number 

appropriate dimensions 

M 1,R 

relative Mach number at sta- 

C P 

specific heat at constant pres- 


tion 1 

F 

sure, ft-lb/(Ib){ 0 R); m/K 

N 

coefficient matrix 

C v 

specific heat at constant volume, 

P 

2 2 

pressure, ib/ft ; N/m 

C 1 

ft -lb/ (lb) (°R); m/K 
coefficient, nondiraensional 

Q 

R 

heat, ft-lb/lb; m 
universal gas constant, ft/°R; 

C 2 

coefficient, sec 


m/K 

C 3 

coefficient, sec 

S 

entropy, ft-lb/(ib)(°R); m/K 

D 

delay matrix 

s 

Laplace transform operator 

d n 

constant coefficients (n = integer). 

T 

column matrix 

E 

£ 

appropriate dimensions 
coefficient matrix 
functional parameter defined in 

t 

V 

time, sec 

velocity, ft/sec; m/sec 


eqs. (100) 

W 

weight flow, ib/sec; kg/ sec 

F 

coefficient matrix 

X 

shock-position column matrix 

G 

shock-position transfer function 

X s 

shock-position, ft; m 


matrix 

X 

space coordinate, ft; m 

S 

gravitational constant, unit 

a 

coefficient, nondimen sional 


force/unit mass 

(3 

coefficient, nondimen sional 

H 

general functional relation 

7 

coefficient, nondimen sional 

h 

enthalpy, ft-lb/lb; m 

A 

perturbation quantity 

i 

V“i 

e 

elements of E matrix 
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V 

functional parameter defined 
in eqs. (100) 

j,k,n 

station number associated 
with variable 

Q 

absolute temperature, °R; K 

r 

reference state 

e 

transportation time, defined 

s 

static condition 


in eqs. (105), sec 

T 

total or stagnation condition 

£ 

functional parameter defined 
in eqs. (100) 

a, &,y 

matrix subscript indicating 
variable location in matrix 

p 

weight density, lb/ft 3 ; 


array 


kg/m 

1,2,3, . . . 

station number associated 

a 

transportation time, defined 
in eqs. (105), sec 


with variable, or variable 
location in matrix array 

T 

transportation time, defined 

Superscripts: 


in eqs. (105), sec 

(-) 

steady- state value of vari- 

U) 

frequency, radian/sec 


able 

Subscripts: 


(~) 

Laplace transform of non- 

av 

steady-state value of Mach 
number averaged over 


dimensional small per- 
turbation variable 


section of subsonic duct 

(■) 

first derivative of variable 

D 

E 

disturbance variable 
exit conditions 


with respect to time 
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APPENDIX B 


SHOCK-POSITION EQUATIONS 

For a stationary normal shock in a steady flow field, the continuity equation across 
the shock may be written as 


P 1 A 1 V 1 = p 2 A 2 v 2 (35) 

where stations 1 and 2 are fixed stations just upstream and just downstream of the nor- 
mal shock as shown in sketch a. 


i / !////) 

j Nt±LLLL 

>/,!////. 

Flow , 

V//7777J 


! 1 

"■Normal shock 


f ////////} 

7777777- 


Fixed stations 1 1 


(a) 


Since shock thickness is quite small, static relations across a stationary shock may 
be applied to a moving normal shock if instantaneous velocities are measured relative to 
the shock and stations 1 and 2 are assumed to be a fixed distance from the shock and, 
hence, moving with the shock. Under these conditions equation (35) becomes 


p l A l (v l ' = p 2 A 2 (v 2 " 


(36) 


When the analysis is restricted to small perturbations in shock position X_, and with 

a 

A 1 ~ A 2’ e( l ua ^ i011 (36) becomes 


*2=^l + Pl7' 1A A 


(37) 


The density relation across the moving normal shock is 
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*1 


6M 


1, R 


m 1r + 5 


where R , the relative Mach number, is defined by 


M 


1, R 




Combining equations (37) and (38) gives 


or 


W t 


s(m? r - l) . 

° w i+pi a i 


M l, R + 6 


(‘♦“"f K R ~ *) ± 


W 2= W 1 


1 + 


a, 


T, 1 


Mi 


M l, R + 5 


(38) 


(39) 


(40) 


(41) 


Linearizing equation (41) and noting that in the steady state X g is zero and Mj R 
(where a bar over a quantity designates the steady -state value of the quantity) result in 


W 


AW. x ' 1 

+ — — 

^ a T M l (l + OJM 2 ^ 2 



(42) 


When stagnation conditions are defined in terms of zero velocity relative to the duct and 
the perfect-gas law is applied at stations 1 and 2, it can be shown that 


P T,2 



(43) 


or 
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P T, 2 " P T, 1] 


P / n 

( 6M 1,R 


J/2 


'© \ 7/2 

T, 2 


6 


M 1,R +5 


0 , 


T,l, 


(44) 


Linearizing equation (44) yields 


AP T, 2 _ AP T, 1 t ? p/ AQ T,2 _ AQ T, 1 


- 35 


( M ! - 0 * 


AM 


1, R 


P T, 2 P T, 1 


© T © T J (nm\ - i)(m^ + 5) M t 


(45) 


where © Tj j = © T| 2 = © T - 

The energy equation across the shock may be written as 


1 2g * 2g 


(46) 


The total enthalpy may be computed from 


hrr = CL©™ = h + 1. 


T " p T 


2g 


(47) 


or 


h, = C„0rr , - — 


‘1 - pT, 1 


2g 


h 2 = C p e T,2 - 


(48) 


Substituting equations (48) into equation (46) 


®m l - ©„ 


gc r 


'T, 1 " V T, 2 " * V 1 " V 2> 


(49) 


Noting that (vj - v g ) = (v 1 - X g ) - (v 2 - X g ) yields 


27 


V 1 ' v 2 = a l M l,R " a 2 M 2, R 


(50) 


but 


_ a lK R + 5 ) 


2 Z,R 6m 


1, R 


( 51 ) 


therefore, 


5a 


V 1 - v 2 = 


6M 


L Kr-‘) 


1, R 


(52) 


Combining equations (52) and (49) gives 


0 = ©m 1 

T, 2 T, 1 


X 5 a 
s 


§ c p 6 m i,r 


L K,h->) 


or 


6, 


T, 2 _ 


X, 


= 1 - 


M i, R - 1 


0rp < 3M-, i / o 

T, 1 1, R T, 1 A +0 .2M] 


yl/2 


(53) 


Linearizing equation (53) and again noting that X - 0 yields 


' &0 T, 2 _ 1\ _ _ 1 


0„ 


Or, 


3a T M 1 


Mj - 1 


/ — 2\1/^ 

(l + 0.2 M ^ J 


AX, 


(54) 


The relative Mach number ^ was defined in equation (39) as 

X 

a l 


which is equivalent to 
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ittti in n i t til ■ i l mu I a I I n j 


l t t 1*9 « III 



X, 


M 


( 1 + 0 - 2M l) 


(55) 


T, 1 


Linearizing equation (55) yields 


(l + 0.2 M^) 1//2 . 

am i,r = am i- , AX S 


(56) 


clri 


At station 1 just upstream of the shock = v i/ a i v j = W^/p-jA^. Thus, 


W, 


M i = 


p l A l a l 


(57) 


Combining equation (57) with the perfect-gas relation and using total pressure yields 


M, 


5a T, 1 W 1 
7gP T, 1 A 1 


3 

(l + 0.2 Mj) = 0 


(58) 


Linearizing equation (58) yields 


AM = 'l( 5 + M l) / AP T, 1 . AA 1 AW I Aa T, 1 1 

1 / <1 


M, 


< M i - ‘) 


P T, 1 A 1 


W 


(59) 


Combining equations (45), (54), (56), and (59) results in 


AP T, 2 _ 6 AP T, 1 7 ( M 1 " ■*■) 

P T 2 7M^ - 1 P T j 7Mj - 1 


. AX g 6(l + 0.2 MjJ AX 


A Lr 




7(m?-i)Aw 1 Aa Tj ^ 


im \ - 1 \ w 


(60) 


where A’/ A is dA/dx at the steady-state shock location divided by the area A. 
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The equation 



can be written as 



Linearizing equation {63) 


ASg 

R 


ASj 

R 


35(Mj 


M l( M l +5 )( 7M " 2 




AM 1,R 


( 61 ) 


(62) 


(63) 


(64) 


Combining equation (64) with equations (56) and (59) results in 


AS 2 AS 1 ?(m^ - l) /AW 1 Aa T? X 7 (m^ - l) AP Tj j 


R 


R 


7M j - 1 


W 


7M 1 ” 1 P T,1 


7 ( M 1 ' X ) 

7M j - 1 


— AX - 
A s 


Mj - 1 


AX, 


M 


/ — 9 \ X /^ a 

i n + 0.2MjJ a T 


(65) 


Equations (65), (60), and (42) relate shock position to perturbations inflow, total pres- 
sure, and entropy. Although valid for both upstream and downstream disturbances, con- 
sideration is hereinafter restricted to downstream perturbations, and these equations 
become 
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( 66 ) 


AP T,2 

P T,2 




A- Lm 


- 0^2 


d^B 
dt L,r, 


AWg 

W 


= c. 


J AX„ 
o s 

dt lZ" 


( 67 ) 


where 


AS 2 _ c L A’ „ c d 

~ - — - ^1^3 — —— 
tt A Uj* at ijrp 


_ i«\ - 1) ] 


c i=— 


7M| - 1 


/ _g\l/2 

„ 5(i + 0.2m2) 

■ — rr^ ' l t 


6a T M t 


1 m| - 1 

c 3 = — — l t 

a T M l (l + 0. 2 uf 2 


( 68 ) 


( 69 ) 


H initial conditions are assumed to be equal to zero, equations (66) to (68) can be written 
in terms of the Laplace transform as 


where 


P 2 - “Cl 


(4X*V) ! 


w 2 = c 3 »x s 


S 2 - C 1 


p = 


ap t (s) 


W - AW < S ) 
W 

g _ AS(s) 
R 


( 70 ) 
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It should be noted that in equations (42), (60), (65) to (68), and (70) the subscripts or 
station numbers refer to stations adjacent to and moving with the normal shock. To use 
the equations in conjunction with other relations describing the subsonic duct, it is desir- 
able to transfer to the primed fixed stations shown in sketch b. 


J/ ////// 

L ... At 1 

'//////// 

V/A 

“1 

^ AL| 

at if/ 

Oft 

Flow 


1 j 

n 

r 

t ! 

\ 

Norma! 
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77777777 

Moving stations 
Fixed stations * 

777 

3 

777777 
■ ' 2 

777 

2 

77777/77 

I 


<b> 

The distance between primed and unprimed stations need only be greater than the 
shock motion. In the steady state, total pressure, flow, and entropy are identical at the 
primed and unprimed stations. Dynamically, they differ only by the time required for a 
disturbance to propagate between the stations, and for small perturbations the total pres- 
sure, flow, and entropy may be regarded as equal at the primed and unprimed stations. 
The subscripts in equations (42), (60), (65) to (68), and (70) can be referred to either 
moving or fixed stations. 
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APPENDIX C 


SUBSONIC-DUCT EQUATIONS 

The relations for the subsonic portion of the inlet can be derived from the continuity, 
momentum, and energy equations. The manipulation required to formulate a usable set 
of equations in terms of the chosen variables (total pressure, flow, and entropy) is out- 
lined in the following section. 

For quasi-one-dimensional flow the continuity equation can be written as 


or 


A (pAv)+ — (pA) = 0 

ax at 


3* + p “ + A |P = 0 

0X 9t St 


From an equation of state of the form 



it can be shown that 


dp = JL dP - JL ds 


and the continuity equation becomes 

3W 3A | Ag SP.Ap 8S,n 

ax at a 2 at c p at 

The 0 dS equations for both total and static pressures are as follows: 

© dS = C_ d© - 

P n 


( 71 ) 


( 72 ) 


( 73 ) 


( 74 ) 


( 75 ) 


( 76 ) 
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(77) 


© T ds = C p d© T 


dP r 


Subtracting equation (76) from equation (77) 


(© T - G)dS = C p d(0 T - ©) 


rdP 


T dP 


P>jr> P 


(78) 


or, noting that © T - © = vv2gC , equation (78) becomes 


dp = -£- dP T -£X dv + ds 


g 


2gC 


(79) 


From the relation W = pAv 


dv = — dW - X dp - ~ dA 
pA p A 


(80) 


Equations (74), (79), and (80) can be combined to yield 


fl - X_ \dp = dP „ - -X- dW + £X_ dA 


Ag 


pv 


Ag 2gC 


dS 


(81) 


from which 9P/3t can be determined and substituted into equation (75). The resulting 
form of the continuity equation is 


_v_ SW / vf _ A SW Agp 9P T + Ap L _ vf_\ 9S _ p 9A = o 


a 


2 at 


ax a 2 p at c 


P V 2a 


2/ at at 


(82) 


For quasi -one -dimensional flow of an inviscid fluid the equation of motion is given 
(see ref. 10) as 



+ P viS,.g® 

ax ax 


(83) 


Multiplying equation (83) by A and summing with the continuity equation (71) multiplied 
by the velocity v give 
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(84) 


v — (pA) + pA — + v A, (pAv) + pAv — = -gA — 

at at sx ax ax 

Noting that pAv is the flow W, equation (84) becomes the momentum equation written in 
terms of flow. 




at 


ax 


ax 


ax 


(85) 


Substituting 3P/3x, determined from equation (79), into equation (85) yields 


3W + v 3W + gAp ^ + PAv? = Q 


at 


ax p T ax 


2C ax 


( 86 ) 


Multiplying equation (86) by (±1 + v/a), adding -a times equation (82), and assuming 
k = 1. 4 give the following pair of equations: 


Apg 


ap 


T L 


ap 


at 


— + (a + v) 


5P t~ | + faw 

axj [a 


+ (a + v) 


Apa ^ 
7R [ 


(M 2 - 2) + M 2 (a + v) 3S 


at 


3Sl - a p M 
ax at 

_ i 


(87) 


Apg 


ap 


T L 


ap rp ap T 

T + <v - a.) T 


at 


ax 




9t 


dx 


= - [(M 2 - 2) + M 2 (v - a) ^§1 - ap — (88) 

7R L St axj at 


The energy equation for an inviscid non- heat -conducting fluid is given in refer- 
ence 10 as 


©/3S + v i§\ = dQ 

\at ax/ dt 


(89) 


or, with no heat addition, 
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( 90 ) 


— + V— = 0 

at ax 

Multiplying equation (90) by 2(M + 1) and noting that v = Ma gives 

(2M + 2) — + 2M(a + v) ~ = 0 

at ax 

Similarly, equation (90) multiplied by 2(1 - M) yields 

(2 - 2M) — + 2M(a - v) — = 0 
3t 3 x 


(91) 


(92) 


Adding equations (87) and (91) and equations ( 88 ) and (92) results in 

_I_ J ^ + (a + v)5 + i[f + (a + v)i 

ap T at ax J Ap [_ at ax _ 

+ aM(M + 2 ) [as +(a + v) as + a 3 A =0 ( 93 ) 

7 R |_ 9 t 3xJ A at 

and 

. a) - A T®* + <v - a) i*l 

axj Ap [at axj 

, a M(M - 2) ras +(v _ a) 9 s] + l 9 A = 0 (94) 

7R [at ax j a at 

Equations (93), (94), and (90) can be used to represent the sub sonic duct. With a change 
in notation this set becomes 
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where 


g 6+P T + J_ £W + aM(M + 2) s^s + a. 3A = Q ' 

ap T at Ap at 7R at a at 


6 

ap T 



J_ £W + aM(2 - M) e~S 

Ap at 7R at 


+ £ M = o 
a at 




9§ + v « 
at as 


= o 


— = — + (v ± a) JL 

at at ax 


( 95 ) 


When slow area change 8A/3t = 0 and small perturbations about an equilibrium condi- 
tion are assumed, equations (95) become 


["— + (v + a) — 1 
|_at axj 

["— + (v - a) — 1 
[at axj 



_ . "| 

+ 7M AW M(M + 2) AS 

M 2 + 5 W M 2 + 5 R _ 

_ 7M AW _ M(2 - M) AS 

M 2 + 5 W M 2 + 5 R 


= 0 


= o 



With the substitutions 


(96) 


(97) 


(98) 


a = 


7M 


M 


av 


av 
4- 5 


J3 - 


M av< 2 + M a») 
M av + S 


and 
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y - 


M av< 2 ~ M av> 
M av + 5 


where M ay is the steady-state Mach number averaged across the section, 
(96) to (98) become 



the set becomes 




w 


R 


AP T _ a AW _ y AS 
P m W R 


V = 


AS 

R 


*- -1 

l+( v + a)Ih=0 

9t axj 

i-+ (V - a) -L~U = oi> 

_at 3xJ I 



v 



>j = 0 


equations 


(99) 


( 100 ) 


(101) 
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For the typical section of subsonic duct illustrated in sketch c, 



the general solution of equations (101) is 


$(M) 

jr(x. 




V 

X - X, 


»?(x,t) = ij^j,t - 


This solution represents a pure delay; thus. 


4(x, s) = exp 


, ['7h <x - x i ) ] i<x i' s) 

d (x, s) = expjljL- (x - Xj)J J &!» s ) 
| (x - Xjjjrjfxj, s) 


t)(x, s) = exp 


H the variable x is set equal to Xg, the length of the segment of duct is L = : 
and equations (103) can be written as 


( 102 ) 


(103) 


2 “ X l> 
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where 


P 2 + G!$ 2 + 0S 2 = e' crs (p 1 + otW 1 + PSj) 

ro . iv \ 

* r«t Te/fv * ^ \ 

P 2 - aw 2 - rs 2 = e TB (Pj - a Wj - rSj) 




(104) 


( 105 ) 
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